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The SU(3) chiral Lagrangian for the lightest octets of mesons and baryons 
is constructed on a spacetime lattice. The lattice spacing acts as an ultraviolet 
momentum cutoff which appears directly in the Lagrangian so chiral symme- 
try remains explicit. As the lattice spacing vanishes, Feynman loop diagrams 
typically become divergent due to inverse powers of the lattice spacing, and 
these divergences get absorbed into counterterms such that the standard re- 
sults of dimensional regularization are obtained. One advantage of lattice 
regularization is that power divergences are seen explicitly. In the present 
work, the octet meson masses, the octet baryon masses and the ttN sigma 
term are all computed explicitly to one loop order. 



I. INTRODUCTION 



Chiral perturbation theory (ChPT) [|IJ is a low momentum effective field theory for QCD. 
For the physics of pions, kaons and eta mesons, ChPT organizes the infinite set of possible 
interactions into a systematic expansion in inverse powers of the chiral scale, A x . Although 
no precise definition of A x is required, it is understood to be 0(1 GeV): 

A x ~ m p ~ ~ 4nF K ~ 1 GeV. (1) 

Powers of 4irF appear as natural suppression factors in the calculation of the loop diagrams 
in ChPT, and the p meson is the lightest hadron which does not appear explicitly in the 
ChPT Lagrangian. The crucial test of ChPT comes from the explicit calculation of observ- 
ables to see whether higher order terms in the expansion really give smaller contributions, 
and the literature contains many such successful examples. [[| 

Since no baryon is light in comparison to A x , it is difficult to include them into ChPT 
without destroying the systematic expansion. Recent work on baryon ChPT is producing 
interesting new suggestions but we will use the traditional solution, heavy baryon ChPT 
(HBChPT) [[|, which is a double expansion in 1/A X and 1/mbaryon- 

In the present work, the chiral Lagrangian is constructed on a spacetime lattice, where 
the inverse lattice spacing plays the role of a momentum cutoff. Because the lattice spacing 
appears directly in the chirally-symmetric Lagrangian, the chiral properties of calculations 
are assured. The lattice-regularized HBChPT Lagrangian is not unique (we will choose an 
isotropic hypercubic lattice) but the continuum limit is unique and is the familiar continuum- 
defined chiral Lagrangian. Since observable quantities cannot depend on the regularization 
and renormalization prescriptions, any valid lattice implementation must reproduce the 
results of dimensional regularization in the limit of vanishing lattice spacing. 

In contrast to dimensional regularization, lattice regularization shows the power diver- 
gences explicitly. These divergences can also be seen with a nonlattice cutoff, but this is 
typically invoked at the time of integration rather than directly in the Lagrangian, so care 
is required to be sure that chiral symmetry is not violated by the cutoff procedure. One ex- 
ample of research that relies upon information about power divergences is found in Ref . || , 
where a nonlattice cutoff was used to discuss the convergence of HBChPT and the scale of 
baryon substructure. 

Lattice QCD simulations are always performed with a nonzero lattice spacing, and they 
typically use unphysical quark masses as well. ChPT is regularly employed for the extrap- 
olations to physical results. Lattice-regularized ChPT offers an analytic way to study 
discretization errors in these chiral extrapolations. Although lattice-regularized ChPT does 
not have a unique Lagrangian, it does allow a determination of the typical size of discretiza- 
tion errors as a function of the quark masses and the lattice spacing. The connection between 
lattice ChPT and lattice QCD has been discussed in Ref. J?J in the limits of strong coupling 
and large N c . 

In section [TT] of the present work, an SU(3) meson Lagrangian is constructed at leading 
and next-to-leading chiral order on a spacetime lattice, and the meson masses and renor- 
malization constants are calculated. This allows a comparison to the work of Shushpanov 
and Smilga || who calculated the quadratically-divergent pieces of F n , and wavefunc- 
tion renormalization in an SU(2) lattice theory. Section [ITT] of the present work contains 
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the meson-baryon Lagrangian. A calculation of the baryon masses and the pion-nucleon 
sigma term is shown to produce the familiar dimensional regularization results as the lattice 
spacing approaches zero. In section fVL the baryon masses and the pion-nucleon sigma term 
are discussed for nonzero lattice spacing, and section [V| contains a concluding discussion. 



II. THE MESON LAGRANGIAN 



When the standard SU(3) ChPT Lagrangian of Gasser and Leutwyler || is written in 
Euclidean spacetime, it takes the following form, 



4? = ^Tr{£ V^V^ - X ] U - X tf}, (3) 

-L 3 £ Tr{ V^VJJVvrfVJJ} + L 4 E Tr{ V ^U}Tr{ X j U + X U ] } 
+L 5 E Tr{ V,U^,U(x ] U + tf X )} - U (Tt^C/ + X U j }] " 



-L 7 (Tr&U - X U ] }) - L B Tr{ X *Ux*U + X U ] X^} 

+iL 9 J2 TV{F* V^V.t/t + F^VjflVJJ) - L 10 J2 Tr{tf F*UF^}, (4) 



where U(x) is a nonlinear representation of the pseudoscalar meson octet, and the current 
quark mass matrix, M, enters via 

X = 2BM. (5) 

External vector and axial vector fields appear within the covariant derivative, V M £/(x), and 
within the field strengths, F^ u (x) and F^(x). 
The chiral transformation of U(x) is 

U(x) -> g(x)U(x)h(x), (6) 

where g(x) G SUr(3) and h(x) G SUl{3). On a spacetime lattice, we introduce the parallel 
transporters, R^(x) G SUr{?>) and L M (x) G SUl(3), which transform as 

R^(x) -> g(x)R^(x)g ] (x + a M ), (7) 
Lp{x) -> h j (x)L II (x)h(x + a M ). (8) 

a M is a Euclidean vector of length a in the // direction, and a is the lattice spacing. The Lie 
algebra- valued vector and axial vector fields, V^{x) and A M (a;), appear in the exponents, 

= exp{-ia£ At (x)}, (9) 
^(x) = exp{-iar M (x)}, (10) 
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where £^(x) = V^(x) — A^(x) and r^ix) = V^ix) + A^ix). In £^ , the following covariant 
derivative is used, 

V«C/(x) = ~ {R p (x)U(x + a,)L\(x) - U(x)} , (11) 
but for all other Lagrangian terms (in , Cy^ > • • • ) a m ore symmetric derivative is used: 
V^U(x) = 1 {^(a;)C/(x + o^a;) - ijt ( a _ - ^(x - o M )} . (12) 

Use of at leading order would produce extraneous poles in the meson propagator 

( "doublers" ) and use of V^ + ) in higher order Lagrangian terms can violate parity. 

The lattice ChPT action is simply obtained by summing the Lagrangian over all space- 
time lattice sites, 

SmP^A^] =a 4 ££ M (x). (13) 

X 

The propagators and vertices required for perturbative calculations can be extracted from 
this action, but if the path integral formalism is used, one must pay particular attention to 
extra meson interactions which get generated by the integration measure, 

8 

DU = e" 5mcasW 11 dvr a , (14) 

a=l 

The definition U(x) = exp{—i\ a 7r a (x)/F} has been employed; A a is a Gell-Mann matrix 
and 7r a (x) is a pseudoscalar meson field. The "effective action" from the measure is [[K]] 

S-M = - - £ Tr In | — 1 , (15) 

<^ = |E^ a (^ ( 16 ) 

b a=l 

where t^ c = —if a b c and f a b c are the structure constants defined by [A a , A 6 ] = 2if abc \ c . These 
measure contributions also exist in the continuum theory, although they happen to vanish 
when dimensional regularization is used. 

Neglecting isospin violation and using mi to denote the up and down quark masses, the 
lowest order pion, kaon and eta two point functions are 



i mm — — yx M 

where 



+ ^X> 2 (?)}. (17) 



x v = J2Bm h (18) 



x K = jB(mi + m s ), (19) 



x v = J^B( mi + 2m s ). (20) 



The meson masses are obtained from the zero of Eq. fllTD , corresponding to the pole in the 
propagator, 

m M = -arcsinh ( aXM \ _ (21) 
a V 2 / 

Notice the existence of a Gell-Mann-Okubo relation, 

Ssinh 2 (^) = 4sinh 2 (^*) - sinh 2 (^) , (22) 

which reproduces the conventional relation as a — > 0. 

At next-to-leading order, the meson masses receive tree-level contributions from C$ and 
from S meas [7r] as well as loop diagrams from the interactions of C^- The only loop topology 
at this order is shown in Fig. [I] and, for example, a charged kaon loop makes the following 
addition to the pion two-point function, 

a -p 1 r /a d 4 p / 2 4 „ . 2 fap v \\ ^ 



3a 2 F 2 J-n/a (2vr) 4 [ K a 2 V V 2 

I + x k) + EX 5 - 4cos(ag M ) - 4cos(ap A1 ) + 3cos(a<^ - ap^)} . (23) 
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Notice that only momenta within the first Brillouin zone can appear in the integral, since 
physics at shorter distances cannot be resolved on the spacetime lattice. The integral is thus 
finite, and only diverges as a — > 0. It is convenient to rewrite the propagator within Eq. ( |23|) 
as the integral of an exponential, using 1/D = / °° dx exp(-xD), which leads to 



6a 4 F 2 



1 + 


a 2 x 2 
2 


W 4 (a 2 x 2 K 


1 + 






K 1 


— -a 2 x 2 K 



(24) 



where 



W n (e 2 ) = | o °°da;/ n (a;)exp|-a;^ + ^| (25) 
and Iq(x) is a Bessel function, 

r* d9 

Iq(x) — J — exp(xcos#). (26) 



By including all of the one-loop diagrams and tree-level pieces, the complete two point 
functions to next-to-leading order are found to be 



MM 



1 ' X M + ^f-E - 2 (^)+( ) E - 2 {a % ) \ (27) 



7 (+) 7 (±) ) 1 a 2 V 2 7 

Z M Z M l a u - ' 
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with 

& = 1 + 24^ + 3^ i l ~ W ^ 

3 



) W 4 (a 2 x 2 K ) - -^L- 2 W 4 (a 2 x 2 v ), (28) 



4 ±} = 1 - + - T2 xlL 5l (29) 

(+) _ 1 , 7 _L 1 (l 3 o 2 \ rxr ./ 2„2' 



^ ( x ~ r 2 ^) + ( x - 2i fl2x2 ) (30) 



4 ±} = 1 - ps(4 + 2a&)L 4 - -^x 2 K L 5 , (31) 

^ ~ i+ 24a 2 F 2 l6F* W * [aXj 

+^r 2 (l - ^M:) ^(a 2 * 2 ,) - ^r 2 W A (a 2 x 2 ), (32) 

4 ±} = 1 " + 2 4)^4 - §- 2 x 2 v L 5 , (33) 

Xt = 4- + 2*k)(L* - 2L 6 ) - A X 4 (L5 _ 2L 8 ) 

+ 2^ + 4^{^^ 
-T^2^(a 2 4) - ^Wk(aV^) - ^W^xt), (34) 
^1 = 4- |^04 + 2 4)(^ - 2L 6 ) - Ax^(L 5 - 2L 8 ) 

+ 24a2F2 + 6a 2 F2^ 4lGX " j 
3 ^ X V 4 (a 2 s*) - ^W 4 (a 2 x 2 K ) - ^W,(a 2 xl), (35) 



64F 2 *y 32j p2 192F 2 

^ 2 = 4 - + 2 ^)(^ - 2^ 6 ) - ^(L 5 - 2L 8 ) + - ^) 2 (3^ 7 + L 

--^ ^(a 2 * 2 ) + { 4^W,(a 2 x 2 K ) + (7 t5 ] ^) 



24a 2 F 2 4a 2 F 2 " v 6a 2 F 2 " v A/ 36a 2 F 2 

_^JL W ^ x l) _ ^L W ^ X 2 ) _ ^Lp^aV). (36) 

The physical meson masses are obtained by using X M instead of xm in Eq. (|2"T|). Notice that 
the pseudoscalar mesons are exactly massless in the chiral limit (mi — m s — 0), indicating 
that the theory does indeed have exact chiral symmetry even for nonzero lattice spacing. 
This feature has been emphasized in the SU(2) case by the authors of Ref. 0. 

The lattice regularized theory gives finite predictions for all observables at nonzero lattice 
spacing. The chiral expansion is in inverse powers of A x . If one chooses to extrapolate to 
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the limit of vanishing lattice spacing, then divergences appear in the loop contributions and 
they must be absorbed into renormalized values of the Lagrangian parameters. With the 
information in Appendix [A], it is a simple matter to show that the continuum limit of the 
lattice theory gives precisely the masses that are familiar from dimensional regularization. 
H Moreover, the logarithmic dependences of the counterterms are found to be 

LJCL/oa) - 2LS(l/o2) - {^(l/n) " 2^(l/ai)} = -3^35 In (^) , (37) 
Ll(l/a 2 ) - 2U s (l/a 2 ) - {L^M) - 2L5(V«i)} = In Q-) , (38) 

3L 7 (l/a 2 ) +L r 8 (l/a 2 ) - {3L 7 (l/ ai ) + L^l/a,)} = -J±— In Q , (39) 

for sufficiently small lattice spacings a% and a 2 . This is precisely the scale dependence that 
is known from Ref. ||, as required, since observables cannot depend on the regularization 
prescription. 



III. THE MESON-BARYON LAGRANGIAN 



The HBChPT Lagrangian is organized as a systematic expansion in the inverse baryon 
mass as well as the inverse chiral scale, A x . This is accomplished by writing the Lagrangian 
in terms of a heavy baryon field, B v (x), instead of the relativistic field, B(x), as follows, 

B v (x) = exp(zm H Bf ■ x )^(l + f)B{x)-, (40) 

where the mass parameter mne is chosen to cancel, or nearly cancel, the octet baryon 
masses. The first few orders in the double expansion of HBChPT are well known [f|], and 
in Euclidean spacetime one finds 

£mb = £mb + Am + £mb + ^mb + higher order, (41) 
£mb = (™o - m HB )Tr (B V B V ) , (42) 

4ffi = E [ Tr (BvV^D^) + X?Tr (B v S,{u^ B v }) + ^Tr (B v S,[u^ B v }) 1 , (43) 

C$ B = ^Tr (B v (v -DvD- D 2 )B v ) - fc^Tr (b v { X+ , B v }) - ^Tr (b v [ X+ , B v \) 

-b Tr(B v B v )Tr( X+ ) + ---, (44) 
4? B = • • • , (45) 

where the omitted terms do not contribute to the present work. m is the leading contri- 
bution to the octet baryon mass that would appear in a relativistic Lagrangian, T> and T 
are the two axial couplings, and = |7s J2v a fiv v v is the Pauli-Lubanski spin vector. The 
matrix B v denotes the baryon octet, 
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B„ 



V 



. J_VO I J_ A 



n, 



2_A 



(46) 



while the pseudoscalar mesons appear within U = £ 2 and x+ — vXs + £x v The hadron 
fields transform under local chiral transformations as 



B v {x) — > o(x)S„(x)o^(x), 
e(x)- ff (x)e(x) t(x) = o(x)e(x)/i(a:). 



(47) 
(48) 



It is convenient to choose v = (0, 0, 0, 1), and to use a nearest-neighbour covariant deriva- 
tive in the time direction, 

aD^B v (x) = B v (x) 



%}(x)R\(x — a 4 )£(:r — a 4 )B v (x — a 4 )^(x — a 4 )i? 4 (a; — a 4 )£(:r) 



— -^(x)L 4 (x — 0,4)%* (x — a4)B v (x — a 4 )£(:r — a 4 )L 4 (:r — a 4 )^(x) 

— (x)R\(x — a 4 )£(:r — a4)B v (x — a 4 )£(:r — a 4 )L 4 (:r — 0,4)^ (x) 

— -^(x)L 4 (a; — a 4 )^(x — a 4 )B v (x — a 4 )^(x — a 4 )_R 4 (x — a 4 )£(:r). 



(49) 



Spatial covariant derivatives appear in C^ B + + . . ., and for these a symmetric definition 
is used (to conserve parity): 

aDjB v (x) = ^(x)Lj(x)£(x + aj)B v (x + a^x + aj)L){x)g{x) 



+ -i\x)Rj{x)i{x + aj )B v (x + aj )Z(x + aj)L](x)g(x) 



+ -S(x)Lj(x)g(x + aj)B v {x + a.j)Z\x + aj )R}(x)Z(x) 

+ \^( x )Rj( x )C( x + aj)B v (x + aj )^(x + a 3 )R){x)i{x) 
1 



£(x)R^(x — aj)ti{x — aj)B v (x — a,j)£*(x — a,j)Rj(x — a,j)£(x) 



nr. — n ■ 



~ g£,( x ) L ]( x ~ a j)i\ x - aj)B v (x - aj)£(x - a,j)Lj{x - aj)(j}(x) 
— -^(x)Rj(x — dj)£,(x — dj)B v {x — dj)i{x — a,j)Lj(x — a,j)^(x) 
~ \^ X ) L )( X ~ a j)£\ x - aj)B v (x - <ij)g(x - aj)Rj(x - aj)£(x). 



One also defines 



(50) 



(51) 



where V^'U(x) is given by Eq. (|12|). 



From £^ B + £^ B , the lowest order baryon two-point function is 



T B b = m HB ~ m o~~ | sm ( a <?4) - 2i sin 2 (j^j } , (52) 
which has a unique zero in the first Brillouin zone, occurring at 

E = —ip4 = - ln{l + a(m — ^kb)}- (53) 
ct 

The physical baryon mass is then whb + E = m + 0(a). Typically the parameter 771hb is 
chosen to equal mo, and then E — at this order in lattice- regularized HBChPT. 

The contribution of £mb to the two-point function is purely tree-level and the contri- 

(3) 

bution at the order of Cl/ B is purely from loop diagrams. The two topologies for loop 
diagrams are shown in Fig. [| The diagram with no internal baryon propagator involves the 
same functions that were used in the previous section for meson masses. For example, the 
contribution of the charged pion loop to the proton two-point function is 

(o) i /•*/« d 4 p J 4 _ . 2 fap v 

x |sin(ag 4 ) — sin(ag 4 — ap 4 ) — 2i sin 2 + 2i sin 2 y ^ 4 - — — ^ j , (54) 

= - 16q 3 F 2 { cos ( Q g4) -ism(aq^)} jl - ia 2 x 2 W A (a 2 xl)^ . (55) 



Because the other diagram in Fig. ^| has an internal baryon propagator, its evaluation is 
somewhat more involved. The calculation is outlined in Appendix |B[ 
As anticipated by Eq. (|53|) , the final results for the baryon masses are 

Trig = TTijjB 4" ~ ln(l + clXb) (56) 
(X 



where 



X B = m - m HB - 2(2x 2 K + x\ )b + 7b^d + ^b T - 



3a 3 F 2 

-Ta^F2 E a BY(a 2 xf) + £ x^ B W 4 (a 2 x 2 ). (57) 

i=-K,K,T) i=TT,K,T] 

Y(e 2 ) = ^(e 2 ) + ^(e 2 ) is defined by Eq. (|A~2|) , and the coefficients and f3 B are given in 
Table |. 

As a — > these baryon masses must be identical to the results of dimensional regular- 
ization. Unlike dimensional regularization, the results of Eq. fl56|) contain cubic and linear 
divergences as a — > 0, but these can be absorbed into renormalized parameters as follows, 
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"* = "*- Jf< je t^P y ^ (58) 

3(P 2 -3:T 2 ) , ^W 4 (0), (60) 

v 128aF 2 v ; 384aF 2 K h K ' 

where a prime denotes differentation with respect to the argument. It is convenient to 
choose 77i jib = ?77q. With reference to Appendix [A|, the a — > limit is easily obtained, and 
is identical to the known dimensional regularization results 11] as required, 



m B -f m - 2(2m^ + mj&o + 7s 6 c + 7^ 2^F 2 ^ 



To conclude this section, consider the pion-nucleon sigma term defined at zero momentum 
transfer via the Feynman-Hellman theorem, 

- dmN (ai\ 
a wN = m -Q^-> ( 63 ) 

where 777 = (m u + m^jl. With lattice regularization and exact isospin symmetry, the 
relation becomes 

a * N = xl {hf) + 1 2W^) XN ' (64) 
leading to the following explicit expression, 

a nN = -2x 2 (2/7 + b v + b r ) 
3a; 2 

+ -^^{W 4 (a 2 xl) + a 2 xlW^a 2 xl) + W 4 (a 2 x 2 K ) + a 2 x 2 K W' A (a 2 x 2 K )} 

+ ^ F2 {W A {a 2 x 2 v ) + a 2 xlW' i {a 2 ^)} 
9x 2 



64aF 2 



*-{V + F) 2 Y\a 2 xl 



x 



2 



bAar z 



Using Appendix the a — ► limit is found to be 

^ - -2m 2 (265 + ^ + &S0 - + - (5X> 2 - + OF 2 ) 

m * m \(V-3F) 2 , (66) 



192ttF 2 



as has been obtained from dimensional regularization fTT|j. 
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IV. THE BARYON MASSES AT NONZERO LATTICE SPACING 



In the previous sections of this work, it has been shown that expressions for the me- 
son masses, the baryon masses and the ttN sigma term are the same in both dimensional 
regularization and lattice regularization in the limit of vanishing lattice spacing. Different 
expressions are obtained when a^O. 

Consider first a lattice spacing that satisfies n/a > A x . Most lattice QCD simulations 
are performed with lattice spacings that satisfy this criterion, but with quark masses that 
are larger than the physical values. If lattice QCD computations are first extrapolated to 
the continuum, then the chiral extrapolations (i.e. the extrapolations of observables from 
the simulated quark masses to the physical quark masses) can use the continuum ChPT 
Lagrangian. Without the initial extrapolation to the continuum, lattice QCD data should 
obey a lattice ChPT Lagrangian instead of the continuum one. In practice, the a^O effects 
of lattice- regularized ChPT should be 0(pa/ir) < 0{p/A x ), where p ~ is a typical 
momentum. This can now be tested explicitly for the observables under discussion. 

For completeness, we also consider a coarser lattice satisfying n/a < A x . In this case, 
Eq. QSBp can be expanded in powers of Xi/(4irF), Xi/m , (ir/a) / (4ffF) and (n/a)/m as 
follows, 

m B = m H B + X B + higher order 

= w!§ + mg + mg + + higher order, (67) 
rrijp = m , (68) 
rri^ = 0, (69) 
mf = -2(24 + 4)bo + Inbv + l T B h T , (70) 

2=7r,K ,7] 

i=-rr,K,ri 

Similar notation will be used for the irN sigma term, 

a nN = afl + a^ ] N + a^ ] N + afl + higher order, (72) 

where 

= *S = 0, (73) 
a { ^ = -2xl(2b + b v + bf), (74) 

and crfjrfj is obtained by subtracting cr^y from Eq. fl65l). Recall that the parameters bo, b-p and 
bp are not dimensionless; they contain an implicit suppression factor due to their position 
in 6 2) 

The four baryon masses plus the ttN sigma term are a set of five observables which 
depend on five parameters: m , bo, bp, bp and V. The other axial coupling is obtained 
from T = gA — "F> ~ 1.267 — T>. To determine discretization errors (i.e. finite-a errors; 
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differences from continuum results), one could hold the parameters fixed and determine the 
a-dependence of each observable, or conversely hold the observables fixed and determine the 
a-dependence of each parameter. In what follows, the experimentally-measured values of 
the five observables will be used to determine the five parameters at each lattice spacing. 

All five observables are linear in mo, bo, bv and bp when the lattice spacing vanishes, 
but the first term on the right-hand side of Eq. (]7ID introduces quadratic dependences for 
nonzero a. This term serves as a reminder that the a / extension of the continuum 
HBChPT Lagrangian is not unique. In fact, the term can be eliminated by adding new 
a-dependent terms to the original (minimal) HBChPT Lagrangian, such as 

A£$ B cx aTr (B V B V ) (Tr X+ ) 2 + . . . , (75) 



with coefficients fixed appropriately. In the present discussion, we will omit the first term 
on the right-hand side of Eq. ([73]) from calculations. The particular discretization errors 
obtained are expected to be representative of a typical lattice ChPT Lagrangian. 

Without the first term on the right-hand side of Eq. (|7l|), T> and JF are easily obtained 
from the Gell-Mann-Okubo relation, 

3 111 
A G mo = ~™k + -m s - -m N - -m H 



(X? 2 - 3^ 2 
64a 3 F 2 



{4Y(a 2 x 2 K ) - Y(a 2 xl) - 3Y(a 2 x 2 v )} . (76) 



With these couplings in hand, bx> can be extracted from mg — m\, and bp from ms — m^- 
Finally, bo is obtained from a W N arid mo from mjv- 

Figure § shows the resulting value of each parameter as a function of the lattice cutoff. 



The "experimental" value of a n N was set to 45 MeV (T^j, and m n was required to satisfy 
Eq. (^). As expected, the renormalized parameters are essentially independent of lattice 
spacing for it /a > A x ; their values in this region are near the dimensional regularized values, 
which are ultimately attained as a — > 0. Significant lattice spacing dependences occur for 
ir/a < 500 MeV. The cubic dependence of the unrenormalized parameter m on the inverse 
lattice spacing is clearly seen. 



V. CONCLUSIONS 

Lattice regularization is a method for introducing an ultraviolet cutoff directly into the 
chiral Lagrangian without destroying the Lagrangian's chiral symmetry. As the lattice spac- 
ing vanishes, lattice regularization represents an alternative to dimensional regularization. 
Unlike dimensional regularization, the lattice theory displays power divergences explicitly. 

In this work, a lattice ChPT Lagrangian for mesons and baryons has been constructed, 
and hadron masses and the pion-nucleon sigma term have been calculated to one-loop order. 

Researchers have had occasion to apply a nonlattice ultraviolet cutoff to ChPT. (See, for 
example, Ref. ||.) However, the spacetime lattice is a particularly convenient implemen- 
tation, due in part to its explicit preservation of chiral symmetry. Another benefit is that 
the lattice spacing appears directly in the Lagrangian, whereas a nonlattice cutoff would 
typically be invoked for each Feynman diagram at the time of loop integration. 
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Lattice regularization is also relevant to lattice QCD simulations, where ChPT is rou- 
tinely used to extrapolate lattice QCD data as a function of quark masses. The present work 
demonstrates how discretization errors arising from the difference between lattice ChPT and 
continuum ChPT can be calculated for a particular lattice ChPT Lagrangian. 
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APPENDIX A: LIMITS INVOLVING W N (e 2 ) AND Y(e 2 ) 

Collected in this appendix are some limits involving the functions, W n (e 2 ) and Y(e 2 ) = 
Y 3 (e 2 ) + Y 4 (e 2 ), defined by 



13 



W n (e 2 ) = J™ dxl^(x) exp j-z + ^ | 



7T (j™^ I e 2 n ' 

n + --$>os0j , (Al) 



- (2tt)» ^ 2 M=1 
F n (e 2 ) =^J o °°dx I%~\x) {I (x) -/£(*)} exp + 



4 - d n I n \ I e 2 n 



cosflj . (A2) 



As e vanishes, the functions of interest approach a finite limit, 

2 



kmW 3 (e 2 ) = W 3 (0) « 0.51, (A3) 
limW^e 2 ) = W 4 (0) « 0.31, (A4) 
limF(e 2 ) = F(0) » 1.43. (A5) 



By making use of the asymptotic behaviour of the Bessel function, 



/ , \ x^oo exp^xj 
Io(x) -> -7=, (A6) 

V 27TX 



the following useful limits are obtained, 



JlMlMLo, (A8) 



0. (A9) 



When comparing lattice regularized results to dimensional regularization, the following re- 
lation is helpful, 



lim <^ v ; ^ I = lim <j — \ + constant. (A10) 



To understand this connection to the logarithm, integrate Eq. ( |A6[ ) to produce the expo- 
nential integral, 



00 



Ei(- X ) = -J^ dt eM t Xt \ (All) 



and then notice that 

z->6 \ ln(x 



Hm m-x^ = L 
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Finally, the calculation of the irN sigma term makes use of the following limits, 

Hb w ^}=~h (A13) 

lim{e^(e 2 )} =0, (A14) 

lim{eF'(e 2 )} = 0, (A15) 

and the identity 

K(e 2 ) = --(2 + -We 2 ). (A16) 
n \ n 



APPENDIX B: A LOOP CALCULATION 

This appendix outlines the calculation of the Feynman diagram in Fig. Q which contains 
an intermediate baryon propagator. For definiteness, the contribution of a charged-pion loop 
to the proton is chosen. 

The diagram represents the following integral, 

Ar(6) • v r /a d *P \ 2 , 4 v- • 2 (o-Pu 



*0j- w / a (27r) 
x |sin(ag 4 + ap^ + it) + 2i sin 
f iV2 



f aq^ + + ie 
V 2 



-i 



[aF (D + F)52S,sm(apJj . (Bl) 

The precise limits of integration deserve some thought. One might consider integrating 
from -min(7r / a,ir / a + q^) through min(7r/a, n/a — g M ) since this would ensure that both the 
meson and nucleon momenta remain within the first Brillouin zone. However, choosing 
v = (0, 0, 0, 1) and working in the proton's rest frame makes q^ suppressed by the inverse 
baryon mass, so the limits of Eq. ( |B"T| ) are accurate to the order we are working. These 
choices also lead to S4 = 0. 

Notice that an "ze" term has been included in the nucleon denominator of Eq. ( |B1~D , so 
that the singularity can be treated correctly. Using 

sin(ag 4 + ie) + 2isin 2 \ — \ = e _e {sin(ag4) + i[e e — cos(ag 4 )]}, (B2) 



V 2 

and a property of the Pauli-Lubanski vector, 



2S , M S , V = ^(v,j,v v - 5^) + [Sp, S v ], (B3) 



one arrives at the following form, 
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X g 1 Sin2( '"' t) { 1+ ( ooshe-ts( a p 4 ) )}- < B4 > 

The piece without sinhe is expressible in terms of Y^a^x 2 ), and the piece containing sinhe 
can be expressed in terms of Y 3 (a 2 xl). The functions Y n (e 2 ) are discussed in Appendix 
The final result is 

ATg = -^r 2 (V + T) 2 {Y 3 (a 2 xD + Y 4 (a 2 xl)} . (B5) 
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TABLE I. Coefficients that appear in the baryon masses. 
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FIG. 1. The only meson loop topology which contributes to a meson mass at one- loop order in 
ChPT. 



(a) 



(b) 



FIG. 2. The two loop topologies which contribute to a baryon mass at leading-loop order in 
HBChPT. 
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(a) 
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FIG. 3. Five parameters, obtained by fitting to the experimental values of gA, &nN and the four 
masses of the octet baryons. The sixth parameter is easily obtained as T = gA — *D. The fit 
to experimental values is redone for each lattice cutoff, it /a. Both the (a) renormalized and (b) 
unrenormalized values of the parameters are plotted. 
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